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ABSTRACT 

An  iterative  technique  for  the  retrieval  of  temperatures 
at  each  of  100  levels  ranging  from  1000  mb  to  0.01  mb  is 
evaluated.   Clear-column  radiance  data  in  the  carbon  dioxide 
channels  of  the  VTPR  of  NOAA  2  are  used  in  inverting  the 
radiative  transfer  equations  to  deduce  the  T(P)  profile.   The 
retrieval  technique  includes  the  computation  of  atmospheric 
transmittances  due  to  three  atmospheric  absorber  masses 
(carbon  dioxide,  water  vapor,  and  ozone)  and  non-homogeneous 
temperature-pressure  effects  along  the  vertical.   The  program 
also  corrects  these  transmittances  for  zenith  path  differences 
between  the  satellite  and  the  retrieval  site  when  the  site  is 
not  directly  below  the  sensor.   A  normalized  wave  number, 
close  to  700  cm   ,  but  decreasing  with  decreasing  values  of 
pressure,  was  used  in  solving  the  Planck  function  for  tempera- 
ture, level-by-level.   The  retrieved  temperatures  were  then 
compared  with  matching  radiosondes  at  15  key  atmospheric 
levels  and  the  resulting  errors  analyzed  on  a  globally  averaged 
and  latitude  band  basis  at  each  of  the  15  levels. 
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I.   INTRODUCTION 

The  quantitative  measurement  of  temperatures  at  known 
pressure  levels  in  the  atmosphere  is  fundamental  to  present 
day  weather  analysis  and  forecasting.   Until  the  launch  of 
the  Nimbus  3  satellite  with  its  Satellite  Infrared  Spectrometer 
(SIRS)  in  April  1969,  free  air  temperature  soundings  were 
obtained  entirely  from  balloons,  aircraft,  and  rockets.   Since 
these  sources  were  readily  available  only  from  populous  land 
areas  and  a  few  ships,  large  areas  of  the  world  were  without 
any  reports.   World-wide  distribution  of  temperature  sounding 
data  was  made  possible  by  the  use  of  satellites.   The  retrieval 
of  a  temperature  profile  by  remote-sensing  satellite  measure- 
ments was  first  suggested  by  Kaplan  (1959) ,  and  has  been  the 
object  of  numerous  research  studies. 

The  Nimbus  3  and  Nimbus  4  satellites  were  equipped  with 
SIRS-A  and  SIRS-B  respectively.   These  first  generation  SIRS 
instruments,  SIRS-A  with  a  vertical  scan,  and  SIRS-B  with  an 
additional  side  scan,  measured  monochromatic  radiances  in  seven 
narrow  regions  or  channels  of  the  15  ym  band  of  carbon  dioxide, 
and  one  in  the  window  region  of  water  vapor  centered  at 
899  cm"  .   Figure  1  (from  Gelman  et  al.  19  72)  shows  the 
weighting  functions  dij/ddnP),  after  Smith  et  al.  (1972), 
of  the  radiometric  emitting  channels  of  the  SIRS-B  instruments. 
A  peak  in  di ./dx  relative  to  the  In  P  axis  depicts  the  layer 
in  the  vertical  contributing  the  maximum  radiance  to  the 
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satellite  channel  sensor.   These  radiances  could  then  be 
related  to  temperature  profiles  through  the  Radiative  Transfer 
Equation  (RTE) 


I±  =  Bj_[T(Ps)]  Ti(Ps)  +   / 

s 


B.[T(P)]  w 


dP 


(1) 


Here,  I.  is  the  spectral  radiance  in  channel  i  (i=l ,2 , . . . , 8) , 
B.[T(P)]  is  the  well-known  Planck  function  for  wave  number  v. 
and  at  temperature  T(P),  t.(P)  is  the  atmospheric  transmittance 
averaged  over  wave  number  i  (where  i  signifies  the  channel 
number  centered  at  v . ) ,  and  P  is  the  pressure-level  for  which 
t- (P)  has  been  computed. 


3 
to 


FIGURE  1.   Derivative  of  transmittance  (t)  with  respect 
to  x(P)  =  P2/7  for  SIRS-B 
(after  Gelman  et  al.  1972) 
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Retrieval  of  temperature  profiles  from  the  SIRS-B  radi- 
ometer was  continued  in  systematic  fashion  by  the  National 
Environmental  Satellite  Service  (NESS)  until  the  radiometer 
output  was  considered  non-usable  on  8  April  1971.   The  number 
of  retrievable  satellite  soundings  collected  during  1969-71 
was  insufficient  to  make  a  marked  impact  on  the  analysis 
problem  over  the  oceans. 

Anticipating  an  increased  demand  for  quantitative  weather 
observations  from  orbiting  satellites,  the  National  Oceanic 
and  Atmospheric  Administration  (NOAA)  in  conjunction  with  the 
National  Aeronautics  and  Space  Administration  (NASA) ,  developed 
a  second  generation  set  of  remote-sensing  radiometers,  improved 
over  the  SIRS-type,  that  would  provide  twice  daily  coverage 
of  most  of  the  globe  and  a  higher  resolution  at  individual 
scan  spots. 

As  a  result  of  these  developments,  the  first  of  the  new 
instruments  was  the  Infrared  Temperature  Profile  Radiometer 
(ITPR)  as  used  on  Nimbus  5.   This  radiometer  carried  four 
carbon  dioxide  channels,  two  window  channels  at  4.3  urn  and 
11.5  urn,  respectively,  and  a  water  vapor  channel  near  20  ym. 
In  1972,  NASA  launched  the  NOAA-2  satellite,  which  carried 
aloft  a  new  remote-sounding  radiometer  for  which  the  ITPR 
had  been  a  prototype.   The  new  instrument,  the  Vertical 
Temperature  Profile  Radiometer  (VTPR) ,  had  been  still  further 
improved  over  the  ITPR  in  its  angular-array  scanning  technique. 
The  VTPR  carried  six  carbon  dioxide  channels,  a  window  channel, 
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and  a  water  vapor  channel.   Observations  made  at  scan  spot 
arrays  in  channels  1,2,.  ..,8  made  it  possible  to  systematically 
retrieve  clear-column  radiances  at  a  higher  resolution  than 
before.   The  new  instrument  was  designed  to  provide  two  tem- 
perature profiles  per  grid  point  per  day  on  a  400  kilometer 
spacing  over  the  oceans,  or  a  global  total  of  6000  temperature 
profiles  per  day.   An  actual  yield  of  72%,  or  about  4000 
profiles  per  day  has  been  realized  [Jastrow  and  Halem  (1973) ] . 
The  use  of  these  VTPR  clear-column  radiances  from  the  NOAA  2 
satellite  and  the  temperature  retrieval  system  described  by 
McMillan  et  al.  (1973)  was  the  basis  for  a  thesis  by  Moran 
(1974) .   The  basic  differences  between  the  Moran  study  and 
previous  work  by  Smith  et  al.  (19  72)  was  the  use,  by  Moran, 
of  standard  atmospheric  carbon  dioxide  transmittances  with  no 
interactive  effects  on  the  channel  1,2,..., 6  transmittances 
from  the  other  atmospheric  constituents.   Moran  also  computed 
radiances  by  the  use  of  a  17-layer  quadrature  scheme  for  inte- 
gration of  Eq.  1  rather  than  by  means  of  the  more  usual  100- 
level  atmosphere.   This  was  done  in  order  to  attempt  to  reduce 
the  trunction  error  involved  with  uncertainties  in  the  integra- 
tion of  the  function  B . [T (P) ] di . /dp  when  a  relatively  crude 
model  for  x.(P)  was  utilized. 

The  first  objective  of  the  present  thesis  was  to  improve 
upon  the  Moran  retrieval  technique  by  incorporating  recently 
updated  state-of-the-art  methods  [Fritz  et  al  (1972) ]  for 
computing  atmospheric  transmittances,  including  the  combined 
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effects  of  water  vapor,  ozone,  and  carbon  dioxide  in  each 
channel,  as  well  as  their  vertical  path  dependencies  on 
pressure  and  temperature.   It  was  also  necessary  to  correct 
these  transmittances  in  each  channel  for  the  applicable  zenith 
angles  of  the  clear-column  spot.   For  this  purpose,  a  special 
data  set  for  March  1973,  including  the  tuned  transmittance 
coefficients  for  the  period,  coupled  with  the  first-guess 
temperature  and  humidity  profiles,  was  acquired  from  NESS. 

A  second  and  equally  important  objective  was  to  test  the 
final  temperature  retrievals  by  comparison  with  radiosondes 
at  island  and  coastal  stations  matching  the  VTPR  scan  spots 
in  both  time  and  spatial  location.   The  retrieved  profiles 
were  compared  at  each  of  the  15  standard  pressure  levels 
(Table  3)  with  the  matching  radiosondes  and  were  analyzed 
both  by  latitude  and  by  pressure  level. 

The  retrieval  system  used  here  has  the  feature  of  inte- 
grating Eq.  1  over  33  atmospheric  layers  at  each  clear-column 
radiance  scan  spot.   This  system  was  incorporated  in  order 
to  minimize  the  truncation  error  in  the  quadrature  scheme 
applied  to  the  function  B . dx . /dp  which  appears  in  the  RTE 
(Eq.  1) .   The  improved  accuracy  in  the  determination  of  the 
temperature  profile  using  this  method  is  one  of  the  main 
thrusts  of  this  thesis. 
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II.   SATELLITE  DATA 

The  NOAA  series  of  satellites  are  modifications  of  the 
earlier  Environmental  Science  Service  Administration  (ESSA) 
satellites.   The  major  change  was  the  elimination  of  all 
vidicon  camera  systems  in  favor  of  radiation-sensing  instru- 
ments, and  in  particular  the  inclusion  of  the  vertical 
temperature  profile  radiometer  (VTPR) . 

A.   NOAA  2  DESCRIPTION 

The  NOAA  2  satellite  orbits  the  earth  at  an  altitude  of 
1464  kilometers  once  during  each  115  minutes.   Figure  2 
illustrates  the  projection  of  seven  orbits  onto  the  earth. 
The  solid  lines  indicate  the  north-to-south  or  "descending" 
portions  of  the  orbits,  which  occur  over  the  sunlit  hemisphere; 
dashed  lines  are  nighttime  "ascending"  portions  of  the  same 
orbits.   Equatorial  crossings  occur  at  0900  (descending)  and 
2100  (ascending)  local  solar  time. 

The  VTPR  instrument  scans  perpendicularly  relative  to  the 
subsatellite  track  from  left  to  right,  when  viewed  downtrack, 
in  23  individual  non-overlapping  steps  covering  ±30.3  degrees 
from  local  nadir.  The  instrument  observes  each  spot  in  eight 
spectral  intervals  or  channels  through  the  use  of  the  automatic 
optical  filter  selector  equipment  of  the  VTPR.  Six  of  the 
channels  are  located  in  the  15  ym  band  of  carbon  dioxide,  one 
in  the  19  ym  region  of  the  rotational  band  of  water  vapor,  and 
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one  in  the  12  pm  atmospheric  window  (see  Table  1) .   After 
the  entire  set  of  8  spectral  measurements  has  been  obtained, 
the  scanning  mirror  rotates  to  the  next  spot. 


Table  1.   VTPR  channel  designators  and  centroid  wave  numbers 
for  filter  set  four  of  NOAA  2  VTPR  instrument 
number  1  (used  in  March  1973) . 


CHANNEL      12345678 
WAVE  NO.   668.5  677.5  695.7  707.2  724.7  747.5  534.5  835.1 
(cm"  ) 

B.   CLEAR-COLUMN  PROCEDURE 

In  order  to  obtain  an  accurate  temperature  retrieval, 
the  "raw"  radiances,  that  is  the  radiances  measured  at  the 
VTPR  instrument,  must  be  adjusted  to  eliminate  the  effects 
of  any  clouds  that  may  be  located  in  the  field  of  view. 
Smith  (19  70)  and  others  introduced  cloud  layer  models  wherein 
the  extent  and  thickness  of  cloud  layers  which  would  give  the 
radiance  as  sensed  at  the  satellite  are  mathematically  intro- 
duced into  the  temperature  profile  calculation.   However,  the 
results  in  the  neighborhood  of  a  single  scan  spot  tend  to 
result  in  non-unique  cloud  heights  and  amounts  which  causes 
uncertainty  in  the  deduced  clear-column  radiances.   The 
method  described  below  was  later  developed  by  Smith  (19  71) 
in  order  to  eliminate  these  uncertainties. 

The  National  Environmental  Satellite  Service  (NESS) 
eliminates  cloud  contamination  through  the  use  of  a  conversion 
program  called  CLRAD,  applied  to  a  large  array  of  contiguous 
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raw  radiance  spots.   Sets  of  eight  scan  lines,  containing 
2  3  scan  spots  are  divided  into  3  subarrays  or  boxes  of  8  by 
8,  8  by  7,  and  8  by  8  scan  spots  as  shown  in  Figure  3.   Ob- 
servations of  raw  VTPR  radiances  are  made  in  each  of  the  eight 
channels  and  at  each  scan  spot.   Figure  4  shows  the  lateral 
scan  geometry  along  each  scan  line,  where  spot  5,  spot  12, 
and  spot  19  in  line  4,  are  the  center  of  box  left,  box  center 
and  box  right,  respectively.   Each  box-center  is  denoted  by 
an  X  in  Figure  3.   At  these  locations  it  is  necessary  to 
determine  the  exact  zenith  angles  to  the  sensor  aboard  the 
satellite  using  the  following  relationship 


SIN(Z  )  =  SIN(N)  2™.  (2) 

n         n        z 


Here  Z   is  the  zenith  angle  at  scan  spot  n  (n  =  5,12,19)  . 

N   is  the  nadir  angle  at  the  same  spots,  a  is  the  constant 
n 

6370  km  (earth's  radius),  and  H  is  the  altitude  of  the 
satellite  (mean  value  of  1460  km) .   The  mean  zenith  angle  at 
spot  12  is  0°,  while  at  spots  5  and  19  it  is  23°47'  left  and 
right  of  center  respectively. 

In  addition  to  the  VTPR  radiances,  measurements  in  a 
Window-channel  Scanning  Radiometer  (SR) ,  sensing  in  10.5  to 
12.5  urn,  permits  determination  of  the  sea  surface  temperatures 
at  each  of  the  scan  spots.   McMillin  et  al.  (p  25-33,  1973) 
outline  a  process  of  extrapolation  of  raw-radiance  to  clear 
column  radiance  within  localized  sub-areas  about  the  scan 
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Figure  3  -  VTPR  scan  pattern  and  data  analysis  array 
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Figure    4.      NOAA-2    VTPR   lateral    scan    geometry   at   raw 
radiance    spots 
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spots.   The  specifics  of  the  clear  radiance  conversion  method 
involves:   (1)  the  specifications  of  a  reasonably  accurate 
value  of  the  first  guess  temperature  profile  T(P)  to  10  mb , 
(2)  an  estimation  of  the  water  vapor  transmittances  in  the 
window  channel  derived  from  the  T(P),  and  (3)  the  determina- 
tion of  the  sea  surface  temperature  derived  from  the  SR.   The 
value  of  the  clear-column  window  channel  radiance  can  then  be 
computed  and  compared  with  the  observed  window  channel  radi- 
ance.  If  the  observed  "raw"  radiance  is  equal  to  or  greater 
than  the  computed  radiance,  the  sounding  is  considered  to 
be  cloud  free.   If  the  computed  radiance  is  greater  than  the 
"raw"  radiance,  the  sounding  is  considered  to  be  cloud 
contaminated  and  the  extrapolation  process  is  necessary. 

To  facilitate  the  extraction  of  the  clear-column  radiance 
set,  "raw"  radiance  values  from  adjacent  spot  pairs  (spots  1 
and  2  with  their  associated  I,  (v)  and  I2  (v)  for  example)  are 

extrapolated  to  I  ,    (v) .   This  can  be  done  if  the  cloud 

c  xear 

amounts  vary  linearly  in  the  field  of  view  and  share  approxi- 
mately the  same  cloud  top  level.   Figure  5  shows  the  patterns 
used  by  NESS  for  pair  comparisons.   Each  8  by  8  box,  so 
analyzed,  will  contain  49  scan  spots,  which  can  be  compared 
with  four  adjacent  spots  for  a  total  possible  yield  of  196 
values  of  clear-column  radiance  I  ,   r(vj)*   Here  i  denotes 
all  channels  other  than  the  window  channel  since  the  window 
channel  gives  rise  only  to  a  possible  49  clear-column  esti- 
mates.  A  proportionately  smaller  number  of  both  radiance 
sets  are  found  in  the  8  by  7  box.   If  a  suitable  number  of 
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the  scan  spots  are  capable  of  giving,  four  adjacent  clear- 
column  values  without  excessive  variances,  all  resulting 
values  are  usable  in  obtaining  the  weighted  clear-column 
radiance  in  the  respective  channels. 


STEP  BY  STEP 
TO    LINE    7 


STEP  BY   STEP 
TO   SPOT  22 


SPOT 


Figure  5.   Scanning  patterns  used  by  the  National  Environmental 
Satellite  Service  (NESS)  for  pair  comparison  in 
determining  VTPR  clear-column  radiance. 
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Research  on  refining  the  accuracy  of  determination  of 
clear-column  radiances  is  actively  being  conducted.   This 
involves  both  instrumental  accuracy  improvements  as  well  as 
the  algorithms  for  the  extrapolations  of  cloud  amounts  to 
zero. 


III.   TRANSMITTANCE  MODELS  IN  THE  VTPR  CHANNELS 

A.   ATMOSPHERIC  TRANSMITTANCE  COMPUTATIONS:   GENERAL 
CONSIDERATIONS 

Atmospheric  transmittances  for  the  six  carbon  dioxide 
spectral  intervals  used  in  the  retrieval  are  fundamental  and 
must  be  computed  accurately.   The  resultant  transmittance  in 
any  given  channel  and  at  any  pressure  level  P (J)  in  the  atmos- 
phere may  be  expressed  as  the  product  of  the  transmittances 
due  to  the  carbon  dioxide  alone  t[C(J)]  and  the  overlapping 
transmittances  due  to  the  integrated  ozone,  OZ (J) ,  and  water 
vapor,  U(J),  from  the  top  of  the  atmosphere  to  the  level  J 
respectively.   Thus 

t(J)  =  t[C(J)]*t[U(J)]*t[OZ(J)]  (3) 

Here  C(J)  denotes  the  total  absorber  mass  of  the  carbon 
dioxide  from  P(l)  =  0.01  mb  to  P (J) .   The  symbol  OZ (J)  de- 
notes the  corresponding  absorber  mass  of  ozone  expressed  in 
NTP  standard  amounts  per  square  centimeter  column  above  level 
J,  and  U(J)  is  the  precipitable  water  vapor  in  the  same  column. 
As  an  example  of  the  transmittances  produced  by  the  program 
discussed  in  this  thesis,  sample  outputs  are  included  following 
the  program  listing. 

Moran  (1974)  used  only  carbon  dioxide  transmittances 
acting  in  the  vertical  column  (i.e.,  t [U (J) ]=t [OZ (J) ] =1 ,  by 
assumption)  and  further  simplified  the  computation  by  selecting 
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a  mid- latitude ,  standard-atmosphere  transmittance  model  for 
carbon  dioxide  in  the  six  listed  channels.   Moran  found  that 
super-adiabatic  temperature  lapse  rates  tended  to  occur 
between  the  1000  mb  and  900  mb  levels.   Consequently  he  found 
it  necessary  to  adjust  the  values  of  the  transmittance  profile 
somewhat  arbitrarily  in  channels  5  and  6,  since  the  trans- 
mittance from  the  surface  in  these  channels  tended  to  be 
somewhat  higher  than  with  the  more  realistic  product 
transmittances  of  Eq.  3. 

Transmittances  in  the  carbon  dioxide  channels  were  based 
upon  the  laboratory  observations  of  the  carbon  dioxide  infra- 
red spectrum  and  resultant  calculations  by  Drayson  (1966) , 
made  on  a  line-by-line  basis  assuming  a  homogeneous  path 
(for  a  fixed  temperature) . 

In  calculations  applied  to  non-homogeneous  vertical  paths, 
Smith  (1969)  ,  calculates  the  carbon  dioxide  transmittances 
corresponding  to  the  U.  S.  Standard  Atmosphere  (1962);  and 
makes  improvements  on  the  usual  Curtis-Godson  modifications 
[Goody  (1964) ]  applied  to  path  transmittances  for  varying 
pressure,  temperature,  and  absorber  masses.   In  addition, 
modified  profiles  corresponding  to  temperature  differences 
from  the  U.S.  Standard  Atmosphere  by  ±10,  ±20,  and  ±30°K 
were  also  computed  by  Smith,  after  Drayson 's  (1966)  laboratory 
calculations.   A  method  for  applying  temperature  corrections 
to  the  absorber  mass  in  each  non-homogeneous  sub-layer  was 
devised  so  that  the  effective  temperature-modified  carbon 
dioxide  transmittances  could  be  adjusted  by  a  variation  of 
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the  Curtis-Godson  method  applied  to  successive  layers,  using 
the  laboratory  data  of  Drayson. 

The  method  utilized  here  for  carbon  dioxide  temperature 
adjustment  along  the  path  will  be  discussed  in  greater  detail 
following  that  for  water  vapor.   This  latter  method  is  some- 
what more  general  due  to  the  variable  mixing  ratio  of  this 
constituent  in  the  vertical  as  well  as  the  greater  complexity 
of  the  individual  line  structure  in  the  water  vapor  infrared 
spectrum. 

B.   WATER  VAPOR  -  MIXING  RATIO  ESTIMATION  PROCEDURE 

In  order  to  determine  the  transmissivity  due  to  water 
vapor  in  each  of  the  VTPR  channels,  the  precipitable  water 
vapor  integrated  from  J  =  01  to  any  arbitrary  level  must  be 
determined.   For  this  purpose  an  estimation  procedure  for 
determining  the  mixing  ratio  ,  W(J) ,  at  each  level  was  used 
to  represent  the  water  vapor  profile.   The  parameter  W(J) 
was  used  because  it  was  found  to  be  determined  from  T(P) 
alone  with  reasonable  skill,  using  the  method  of  Weinreb 
and  Crosby  (1973) ;  and  also  because  it  enters  directly  into 
the  calculation  of  the  water  vapor  transmittances  in  the 
eight  VTPR  channels. 

The  calculation  of  the  mixing  ratio  W(J)  may  be  derived 
from 

W(J)  =  SAT(T,J) *H(J)  (4) 


28 


H(J)  is  the  relative  humidity  and  is  an  unknown  function  of 
pressure.   The  use  of  the  v/ater  vapor  channel  (v  =  535  cm   ) 
was  not  attempted  for  retrieval  in  this  thesis.   Furthermore, 
SAT(T ,J) ,  the  saturated  mixing  ratio  is  a  function  of 
temperature  T(J)  and  pressure  P(J)  only,  and  is  calculated 
by  the  GOFF-GRATCH  formula  [List   1963] 


log  e   =  -7.90298 (T_/T  -  1)  +  5.02808  log(T/T)        (5) 

-  i.asiexio-^io11-34^1"1/^'  -  i] 

+  8.1328xl0-3[10-3-49149(T^T-1>  -  1]  +  log  e 

^   ws 


Here  e   is  the  saturation  vapor  pressure  over  a  plane  water 
surface , 

T   is  the  temperature  in  °K, 

T   is  the  steam  point  temperature  373.16  °K,  and 

e    is  the  saturation  vapor  pressure  of  water  at  the  steam 
point  temperature. 

The  value  of  SAT(T,J)  is  then  computed  using  the  relationship 


e 

SAT(T,J)  =  0.622(p/T,W — )  (6) 

P(J)  -  ew 


Since  no  equation  directly  relates  the  relative  humidity 
to  the  temperature  profile,  the  mixing  ratio  was  estimated 
from  the  saturated  mixing  ratios  at  a  specified  set  of  eleven 
levels  (Table  2)  by  a  least-squares  formulation  based  solely 
upon  these  eleven  saturated  mixing  ratio  predictors  after 
Weinreb  and  Crosby  (19  73) .   The  method  formulates  the  mixing 
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ratio  W(J)  at  the  40  levels  J  =  61 , 62  ,  . . .  ,  100  and  is  given 
in  the  form 


11 
W(J)  =  W(J)  +  I      CREG(I,J) [SAT(I)  -  T(I)]  (7) 

1=1 


where  W(J)      is  the  sample  mean  mixing  ratio  at  level  J  over 

1100  sample  soundings  that  were  used  to  determine 
the  regression  coefficients, 

CREG(I,J)  is  a  set  of  predetermined  regression  coefficients 
calculated  for  each  J-level  from  the  1100  sample 
soundings , 

SAT (I)    is  the  saturation  mixing  ratio  at  the  level  I, 
where  I  =  1,2,...  ,11  is  a  subset  of  the  40 
J-levels  that  is  found  to  give  the  best  results 
(see  Table  2) ,  and 

T(I)      is  the  mean  saturation  mixing  ratio  at  level  I 
from  the  1100  test  soundings. 


Table  2.   Eleven  levels  used  for  predictor  levels  in 
water  vapor  retrieval 


P(J)  mb 

209.94 
245.37 
299.01 
412.26 
509.93 
699.03 
753.17 
810.25 
870.35 
933.57 
1000.00 


Thus  the  method  uses  an  eleven  level  saturated  mixing 
ratio  profile,  SAT (I),  and  converts  it  into  a  40-level  mixing 
ratio  profile,  W(J),  J  =  61 ,  62 ,  .  .  .  ,100  .   The  total  100-level 
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profile  is  then  converted  from  the  40-level  by  using  the  cubic 
power-law  profile  due  to  Smith  (1966) 

W(J)  =  W(61)  t|^]3  (8) 

J<61 

The  mixing  ratio  is  then  used  in  determining  the  water  vapor 
transmissivity . 

C.   ADJUSTMENT  OF  WATER  VAPOR  MASS  IN  A  NON-HOMOGENEOUS 
ATMOSPHERE 

Assuming  that  the  transmittance  by  water  vapor,  t[U(J)], 

in  a  spectral  interval  of  the  VTPR  is  a  known  function  of  the 

quantity  of  the  precipitable  water  vapor  U,  temperature  T, 

and  pressure  P,  an  accurate  approximation  to  t[U(J)]  in  the 

form  of  a  least  squares  fit  polynomial  can  be  calculated. 

Weinreb  and  Neuendorfer  (1973)  proposed  such  a  best-fit  model 

which  in  effect  becomes  an  analytic  function  for  calculating 

homogeneous  path  transmittances.   In  the  Weinreb-Neuendorfer 

method,  an  inhomogeneous  atmosphere  is  treated  as  a  sequence 

of  homogeneous  layers,  each  having  a  constant  temperature, 

pressure,  and  absorber  mass.   The  transmittance  model  provides 

for  a  homogeneous  path  as  a  known  function  of  T,  P,  and  U, 

where  U  is  the  effective  water  vapor  mass.   If  t(J)  is  the 

transmittance  in  the  spectral  interval  of  the  instrument  at 

level  J,  then 

t(J)  =  t[P(J)  ,T(J)  ,U(J)]  .  (9) 
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The  procedure  involved  in  defining  t  [  (J) ]  at  each  level 
as  depicted  in  Fig.  6,  uses  the  functional  expression,  Eq.  15 , 
for  layer  1  without  modifying  the  precipitable  water  vapor 
mass,  U,  from  its  initially  specified  value  based  upon  the 
homogeneous  values  P(l),  T(l),  and  U(l).   This  gives  rise  to 
the  following  relationship  for  layer  1 


T(l)  =   [P(D  ,T(1)  ,U(1)] 


(10) 


U(l)  E  U(l) 


LAYER 
1 
2 


LEVEL 


Figure  6 . 


n-1 


I"   (SURFACE) 

Multi-level  atmosphere  showing  P(J),  T(J),  and 
U(J)  structures. 


In  proceeding  to  level  2 ,  the  transmittance  through  the  second 
layer  is  affected  by  the  new  homogeneous  characteristics  P(2), 
T(2) ,  and  its  water  vapor  increment  AU2 ,  and  is  no  longer  along 
the  x [P(D ,T(1) ,U(1) ]  transmittance  curve  defined  by  the 


32 


conditions  in  layer  1  (points  A,-B,  in  Fig.  7)  .   It  is 
necessary  to  determine  the  adjusted  value  of  water  vapor  mass, 
denoted  V  in  Fig.  7 ,  with  which  the  radiant  transmittance 
enters  the  new  homogeneous  layer  2.   This  is  accomplished  by 
backing  up  on  curve  2  to  U  =  V  by  solving 

t(1)  =  x[P(2)  ,T(2)  ,V]  (11) 

for  V.   The  solution  for  V  is  accomplished  by  the  Newton- 
Raphson  iteration. 

This  corrected  starting  water  vapor  mass,  V,  with  which 
to  enter  layer  2,  has  been  deduced  using  curve  2  in  Fig.  7. 
For  the  purpose  of  computing  the  transmittance  at  level  2 ,  the 
new  total  water  vapor  mass  is  then  taken  as 


W  =  V  +  AU2  (12) 


where  AU~  is  the  known  increase  in  precipitable  water  vapor  in 
layer  2.   This  moves  the  transmittance  value  t  [P  (2) ,T  (2)  ,W] 
down  to  point  B2  on  curve  2  where  the  value  of  t(2)  is  then 
taken  as 

T(2)  =  t[P(2)  ,T(2)  ,U(2)  =W]  (13) 

In  proceeding  downward  to  level  3  through  the  next  homogeneous 
layer  of  characteristics  P(3),  T(3),  AU_,  the  corrected  amount 
of  water  vapor  at  level  2  is  revised  from  U(2)  to  V  which 
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Figure  7.   Transmittance  curves  for  atmospheric  layers  1  and  2 
illustrating  modification  to  effective  U(J)  in  a 
non-homogeneous  atmosphere. 
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becomes  the  level  2  starting  point  of  water  vapor  mass. 
Hence  the  new  function 

x(3)  =  t[P(3)  ,T(3)  ,U(3)  ] 

(14) 
U(3)  =  V  +  AU3 

is  solved  for  the  new  starting  mass  at  the  top  of  level  3. 
The  procedure  continues  as  a  series  of  Newton-Raphson  itera- 
tions to  correct  the  effective  water  vapor  masses  for  the 
radiant  transmittance  upon  succeeding  layers.   This  procedure 
is  followed  through  each  level  J  =  2, 3,... ,99  bounding  the 
top  of  the  9  8  layers  below  J  =  2. 

The  procedure  is  shown  by  Weinreb  and  Neuendorfer  to  be 
exact  for  a  Beer's-Law  monochronatic  transmittance  function. 

D.   WATER  VAPOR  TRANSMITTANCE  MODEL 

The  water  vapor  transmittance  calculations  use  the  method 
of  Weinreb  and  Neuendorfer  (1973)  applied  to  a  modification 
of  the  polynomial  representation  of  water-vapor  transmission 
of  Smith  (1969) .   The  basic  equation  for  the  line  transmit- 
tance of  water  vapor  implies  that 


14 

Z 

1=1 


ln{-ln  t  [U(J) ]}  =  Z        C1X1  (15) 
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Here 


X1       1  Xg  -  x4x? 

x2  =  (o.i)in~~      x9  =  x3x4 

X3  =  ln(P/1000)       X1Q  =  X2X? 

X4  =  ln(T/273)        X^  =  X4X? 

2 
X5  "  X2X3  X12  =  (X4) 


X6  "  X2X4  X13  ~  X3X6 

2 

X~    —     \X.y)  ~\    A  "3    7 


and  CT  is  a  set  of  coefficients  that  were  predetermined  by 
fitting  Eq.  15  by  the  least-squares  method  using  line-by-line 
calculations  of  the  transmittances  for  130  homogeneous  (con- 
stant pressure  and  temperature)  paths.   From  Eq.  15  it  follows 
that 

14  MCI 

t[U(J)]  =  exp[-exp(  Z  CTXT)  ]  ^D; 

1=1  1  -1 

The  values  of  X-,  through  X,  .    are  directly  calculated  only 
if  U  is  known.   However,  X2  must  be  calculated  indirectly 
since  U,  the  water  vapor  absorber  mass,  is  essentially  a  sum 
of  scaled  values  of  AU  between  successive  J-levels.   Newton's 
method  was  used  to  solve  for  the  new  value  of  X2 ,  and  subse- 
quently other  X.  =  X.X2,  using  the  first  iteration  of 
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uj    +  AU  =  Up    as  the  first  guess  to  arrive  at  the  effective 

(2) 

U2    at  level  2,  etc.  on  to  the  subsequent  J-levels. 

Water  vapor  transmittances  in  the  infrared  region  have 
been  found  to  contain  a  spectral  feature  not  previously  ex- 
plainable by  the  classical  line-spectral  approaches.   Bignell 
(1970)  found  a  continuum  absorption  from  the  edge  of  the  20 
ym  rotational  band  which  overlaps  the  lines  in  the  water-vapor 
and  carbon-dioxide  absorption  regions .   The  continuum  absorp- 
tion coefficient  was  found  to  be  due  to  the  presence  of  the 
water  vapor  dimer  (H20)2.   The  dimer  absorption  coefficient  in 

the  various  spectral  regions  was  found  to  be  proportional  to 

-5 
[T(J)]    (see  Eqs.  17  below).   Thus,  it  is  necessary  to  add 

the  dimer  correction  to  Eq.  15  to  fully  account  for  the  com- 
bined transmittances  of  water  vapor  in  all  channels.   The 
dimer  contribution  at  level  J  to  the  ln[-lni(J)]  of  Eq.  15 
due  to  water  vapor  is  given  by  the  adjusted  dimer  mass  profile 
in  terms  of 


DIMER  (J)  =  DIMER(J-l)  +  tA(J-l)+A(J)  ]  [P ( J) -P (J-l)  ] 


A(J)  =  P(J)  [W(J)]2  [30jj^] 


(17) 
5 


The  actual  transmittance  due  to  the  effective  dimer  mass 
in  each  channel  is  then  given  by 

In [-In  x (DIMER)]  =  B (I) *DIMER(J)  (18) 
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where  B(I),  I  -  1,2,. ..,8  is  a  set  of  8  constant  absorption 
coefficients,  the  values  of  which  have  been  given  for  each 
channel  after  Bignell  (1970).   The  combined  water-vapor 
transmittance  in  each  channel  then  becomes 

T[U(J)]  =  exp{-exp[DIMER(J)*B(I) ])  t(J)  (19) 

where  t  (J)  is  the  value  computed  using  Eq.  16. 

The  zenith  angle  correction  for  water  vapor  transmittance 
is  included  in  DIMER(J)  and  U(J)  of  Eqs.  9  and  17  by  making 
the  following  path-corrections  for  mass: 

DIMER(J)  =  DIMER(J) *PATH, 

U(J)  =  U(J)*PATH,  (20) 

PATH  =  sec  Z . 

Here  sec  Z  is  the  secant  of  the  zenith  angle  from  the  point 
representing  the  center  of  the  box  to  the  radiometer  (see 
Figure  4) . 

E.   CARBON  DIOXIDE  TRANSMITTANCE  MODEL 

A  polynomial  expression  for  the  carbon  dioxide  trans- 
mittance t[C(J)]  in  terms  of  logarithmic  functions  similar 
to  Eq.  15  has  been  shown  to  be  feasible  by  Smith  (1969) . 
However,  the  absorber  mass  of  carbon  dioxide  is  proportional 
to  pressure,  and  little  dependence  upon  higher  powers  of  the 
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ln(P/1000)  term  has  been  found  necessary  to  fit  line-by-line 
computations  of  transmissivities  with  increasing  pressure  in 
non-homogeneous  atmospheres. 

However,  the  carbon  dioxide  polynomial  expansion,  similar 
to  Eq.  15  for  the  water  vapor  model,  has  been  found  to  be 
subject  to  path  correction  primarily  for  the  temperature  varia- 
tion along  the  vertical.   The  method  of  Weinreb  and  Neuendorf er 
(1973)  was  shown  by  Fleming  (1974)  to  be  converted  into  the 
following  form 


t[C(J)]  =  t[C(J-1)  ]  [a(J)+B(J)AT(J)+y(J)AT2(J)  ] 


+  6i{T[C(J-l)]-T*(J-l)} 


(21) 


where  x[C(J)]  is  the  transmittance  due  to  carbon  dioxide  at 
level  J: 

AT (J)  is  the  temperature  difference  between  the  U.  S. 
Standard  Atmosphere  (1962)  and  the  average 
temperature  over  the  layer  bounded  by  the 
pressure  levels  J  and  J-l  (J  =  2 , 3,  .  . . , 100)  , 

a,B/Y  are  a  set  of  3  by  100  transmittance  factors  that 
are  tuned  specifically  to  the  filters  of  the 
instrument  in  use  for  the  data  period, 

0.  is  a  set  of  six  adjustment  factors  applicable  to 

systematic  corrections  for  temperature  varia- 
tions from  the  initial  standard,  and 

T*(J)  is  a  correction  term  approximated  after  the 
Newton-Raphson  method  from 

2 

*,_.  _  A(J)+B(J)  AT  (J)    D(J)  [AT(J)  ]  ,99* 

T*(J) Am/  T^1  V +   T~  ^  ' 

AT(J  X)  [AT(J-1)T 
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where  for  the  initial  level  J  =  1,  A,  B,  D  are  defined  as 

A(l)  =  a(l) 

B(l)  =  Bd)AT(l)  (23) 

2 
Dd)  =  tg   a(if  +Y(1)][AT(1)] 

and  for  subsequent  layers,  J  >  1,  the  values  of  A,  B,  and  D 
are  updated  iteratively  by 


A(1)  (J)  =  A(2)  (J-l) 


B(1)fJ>  -  B(2)  (J-l)AT(J) 

B    (J) AT(J-l) (24) 

D(1)  (j)  =  D(2)  (J-D  [AT(J)]2 

[AT(J-l)]2 

D(2)(J)  ,  D(1)(J)[AT(J)]2a(J)  +a(1)(j)y(j)[at(j)32 

[AT (J-l) ]^ 


B(2)(J)  =B(1)(J)a(J)  +  A(1)  (J)3(J)AT(J) 


A(2)  (J)  =  A(1)  (J)a(J) 


all  of  which  must  be  computed  in  the  given  order. 

F.   OZONE  TRANSMITTANCES 

The  ozone  transmittances  x [OZ ( J) ]  are  minor  contributions 
appearing  in  thie  form  of  Eq.  3  for  each  carbon  dioxide  channel 
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They  are  constant  in  each  channel  at  a  given  level,  varying, 
however  from  level  to  level.   Thus  the  effect  of  t  [OZ  (J) ] 
is  expressible  in  the  form  of  a  6  by  100  array,  one  member 
of  the  array  for  each  carbon  dioxide  channel  at  each  J-level. 
The  channel  values  of  t  [OZ (J) ]  used  here  are  listed  in  Appendix 
B.,  as  extracted  from  the  NESS  CLRAD  program. 

The  t  [OZ  (J) ]  values  in  Appendix  B  have  not  been  path 
corrected  and  are  based  upon  globally  averaged  NTP  absorber 
paths  of  ozone  in  the  carbon  dioxide  channels. 

G.   TRANSMISSIVITY  ADJUSTMENT  FOR  ZENITH  ANGLE 

When  the  transmittance  in  the  spectral  interval  of  the 
VTPR  has  been  calculated,  an  adjustment  must  be  made  to 
account  for  the  path  differences  that  occur  with  various 
zenith  angles.   The  VTPR  is  designed  to  measure  clear-column 
radiances  from  points  which  subtend  a  zenith  angle  of  23°47' 
left,  0°,  or  23°47'  right  of  the  satellite  track.   These 
angles  correspond  to  the  center  of  the  three  boxes  shown  in 
Figure  3.   The  adjustments  to  the  calculated  transmittances 
must  be  made  for  these  basic  angles  and  any  small  deviations 
from  them. 

The  principle  of  the  correction  stems  from  the  fact  that 
if  the  zenity  angle  is  other  than  zero  degrees,  the  atmospheric 
absorbers  will  have  more  effect  due  to  the  geometry  involved 
(Figure  4)  with  the  longer  path,  and,  therefore,  the  trans- 
missivity  in  the  layer  will  be  decreased  level  for  level. 
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The  correction  for  both  the  carbon  dioxide  and  ozone 
transmittance  at  level  J  is  taken  from  a  section  of  the  NESS 
CLRAD  program  [McMillin (1974) ]  and  is  in  the  form  of 


T(J)  =  t(J)A  (25) 


A  =  (PATH)B(J)  (26) 


where  PATH  is  the  secant  of  the  zenith  angle,  and  B(J)  is  a 
polynomial  function  for  carbon  dioxide  alone.   For  the  ozone 
path  corrected  transmittances  B(J)  =  1  for  all  atmospheric 
levels.   B(J)  is  a  function  of  ln[P(J)]  with  terms  in  ln[T(J)] 
up  to  quadratic  (and  a  product  of  these  logarithmic  functions) , 
with  a  set  of  five  predetermined  constants  for  each  channel. 
These  constants  have  been  tuned  for  the  particular  VTPR 
instrument  in  use  during  the  period. 
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IV.   RETRIEVAL  OF  TEMPERATURE  PROFILES 

A.   INPUT  DATA  USED  FOR  TEMPERATURE  RETRIEVAL 
1 .   First  Guess  Temperatures 

In  addition  to  the  climatological  soundings,  T(P), 
from  the  U.  S.  Standard  Atmosphere  Supplements  (19  66)  (see 
Appendix  B) ,  the  retrieval  program  requires  a  first  guess 
temperature  profile  in  the  box  center  for  each  VTPR  clear- 
column  radiance  set.   The  first  guess  temperatures  are  of 
two  basic  types  depending  upon  the  latitude  of  the  sounding. 
For  tropical  latitudes,  from  the  equator  to  18  degrees  North, 
radiosonde  data  is  averaged  over  all  latitudes  and  longitudes 
to  provide  a  single  climatological  first  guess  profile.   For 
the  other  northern  hemisphere  retrievals,  the  first  guess  pro- 
file comes  from  an  application  of  the  National  Meteorological 
Center  (NMC)  prognostic  model  for  a  period  of  12  to,  at  most, 
18  hours  from  the  preceding  complete  analysis.   The  profile 
is  extrapolated  in  time  and  space  to  the  VTPR  clear-column 
sounding  location. 

The  first  guess  temperatures,  as  read  into  the  program, 
give  temperatures  at  the  15  mandatory  levels  shown  by  an  aster- 
isk (*)  in  Table  3,  from  100  mb  to  10  rob.   The  U.  S.  Standard 
Atmosphere  Supplement  climatology  value  is  adjusted  for  latitude 
and  month  and  then  fitted  to  the  first  guess  profile  at  the  10 
mb  level.   The  adjustment,  defined  as  first  guess  (10  mb)  minus 
standard  (10  mb) ,  is  then  added  to  all  values  of  the  standard 
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Table  3.   56  level  atmosphere  with  merged  climatology  and 

15  level  (*)  first  guess  field  for  sample  sounding 

LEVEL  PRESSURE  (mb)  TEMPERATURE  (°C) 

1  0.01  -87.130 


23  1.00  -04.311 


34  10.00*  -42.300 

35  15.00  -45.300 

36  20.00*  -49.925 

37  25.00  -53.175 

38  30.00*  -56.800 

39  40.00  -5S.900 

40  50.00*  -63.000 

41  60.00  -65.600 

42  70.00*  -68.200 

43  85.00  -69.250 

44  90.00  -69.600 

45  100.00*  -70.300 

46  150.00*  -63.400 

47  200.00*  -58.200 

48  250.00*  -48.400 

49  300.00*  -38.900 

50  400.00*  -23.600 

51  500.00*  -11.800 

52  600.00  -03.500 

53  700.00*  4.800 

54  850.00*  13.600 

55  900.00  16.267 

56  1000.00*  21.600 
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atmosphere,  above  10  mb,  and  the  profile  is  expanded  to  the 
56  levels  shown  in  Table  3. 

The  56  level  temperature  profile  is  then  expanded  to 
the  100  J-levels,  where  the  pressure  at  level  J  is  given  by 


7/2 
P(J)  =  0.01[1+ (J-l)  (0.26087836) ]         (27) 


where  J  =  1,2,... ,100  and  defines  the  pressure  from  0.01  mb 
to  1000  mb. 

2 .   Other  Input  Data 

The  remaining  input  data  required,  having  given  the 
VTPR  first  guess  sounding,  includes  the  zenith  angle,  sea 
surface  temperature,  and  clear-column  radiances.   This  data, 
along  with  all  of  the  verifying  radiosondes  and  first  guess 
temperature  profiles  to  10  mb  was  kindly  provided  by  NESS  in 
the  form  of  an  archival  tape. 

B.  MATHEMATICAL  DEVELOPMENT  OF  THE  TEMPERATURE  RETRIEVAL 

Given  a  non-scattering  atmosphere  in  thermodynamic  equili- 
brium, the  clear-column  spectral  radiance  at  the  top  of  the 
atmosphere  in  each  channel  is  related  to  the  vertical  tempera- 
ture profile  and  the  absorber  mass  structure  by  the  RTE , 
Eq.  1. 

The  right  hand  side  of  Eq.  1  is  composed  of  an  earth's 
surface  radiation  term,  and  an  atmospheric  contribution  which 
is  integrated  from  the  surface  of  the  earth  (1000  mb)  to  the 
top  of  the  atmosphere  (0.01  mb) . 
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When  the  100  J-levels  developed  in  Eq.  27  are  used,  Eq.  1 
can  be  rewritten 

J=1  dT  (J) 

I.  =  B.  [T(100)  ]t-  (100)  +   /     (B.[T(J)]    i^.,_ 

1   x       1  j=ioo  x  ~dj — )dJ 

(28) 

The  Planck  radiance  function  used  in  Eq.  2  8  at  level  J  and 
in  channel  i  is 

(Cl)*v,3 

VT(J)]  "  C2^T <29> 

exP[Tljf]  "  1 

where  V.  is  the  wave  number  in  channel  i, 

-5       2     -1      -1 
Cl  is  the  constant  1.9061x10    erg  cm  sec    ster   ,  and 

C2  is  the  constant  1.43868  cm  °K. 

In  order  to  reduce  the  effect  of  truncation  error  in  the 

integration  function  of  the  RTE,  and  increase  overall  stability, 

the  100-level  atmosphere  was  reduced  to  a  33-layer  atmosphere 

as  shown  in  Fig.  8.   In  integrating  the  RTE  over  the  100  level 

atmosphere,  Simpson's  rule  was  applied  over  each  sublayer  K 

centered  at  the  fractional  J  levels  of  the  K-layered  atmosphere 

(Fig.  8) .   The  main  virtue  of  the  Simpson  quadrature  scheme 

using  the  weighting  factors  of  1,  4,  1  over  each  layer  (e.g., 

J  =  01,  02.5,  04),  lies  in  the  small  discretization  error, 

5 
(1.5)  /90  times  the  fourth  derivative  with  respect  to  J  of  the 

integrand  function  of  Eq.  28,  in  substitution  of  the  finite 

summation  for  the  Riemann  integration.   The  determination  of  T 
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Figure  8.   Thirty  three  K-layer  atmosphere  model  used  for 

the  integration  of  the  Radiative  Transfer  Equation 
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at  levels  J  =  02.5  etc.,  is  accomplished  by  interpolation  of 

2/7 
the  T(P)  using  the  P     pressure  levels,  defined  by  Eq.  27. 

Therefore,  Eq.  2  8  can  be  rewritten 


33 
I.  =  B.  [T(100)]T.(100)  +  I      B.  (K)At.(K)  ,  (30) 

K=l 


where 


Bi(K)  =  |{Bi[T(JK-1.5)  +  4Bi[T(K)]  +  B± [T (JR+1 . 5) ] } 


(31) 


is  the  layer-mean  Planck  function  for  layer  K  employing 
Simpson's  rule.   Moreover  At. (K)  now  spans  three  J  levels 
(Fig.  8) , 


At±(K)  =  Ti(JK-1.5)  -  xi(JK+1.5)  (32) 


where  J„  is  the  value  of  J  at  the  center  of  layer  K. 

Using  the  inversion  method  of  Smith  (1970)  as  modified  by 
Martin  (1974) ,  the  retrieval  of  mean  layer  temperatures  is 
determined  from  clear-column  radiance  values  as  follows: 

1.  Determine  the  first  guess  radiance  value  in 
channel  i,  I.   ,  using  Eq.  30. 

2.  Form  the  residual  by  comparing  the  clear-column 
radiance  with  the  computed  value  from  Eq.  30,  using 


"I1'  "  ^-CLEAR  "  ij1'  (33) 
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3.   If  the  residual  is  not  close  to  zero  the  first 

guess  temperature  profile  T    (p7  )  must  be  improved  to 

(2)  K 

T    (PT  )  in  accordance  with 

JK 

B.  [T(2)  (P_  )  ]  =  B.  [T(1)  (PT  )  ]  +  R.(1)  (34) 

1        JK       X        JK  1 


This  tends  to  move  the  temperature  profile  closer  to  the  final 
desired  convergent  profile.   If  n  steps  are  required  to 
converge  the  residual  close  to  zero,  the  final  step  is 


B.  [T(n+1)  (P   )  ]  =  B.  [T(n)  (P_  )  ]  +  R.(n)  (35) 

1  JK       1        JK       1 


where  finally 


R.(n) 

-= <  10  4  (36) 


Ii-CLEAR 


When  a  convergent  Planck  function  is  determined  for  each 
channel,  i,  and  at  each  level,  J,  the  temperature  at  that 
level  can  be  determined  by  solving  Eq.  29  for  T(J).   The 
value  of  the  Planck  function  in  this  case  must  be  a  weighted- 
mean  Planck  using  the  six  values  of  B.(J),  i  =  1,2,.  ..,6,    at 

each  J-level.   Figure  9  shows  the  values  of  transmittance  in 

2/7 
each  of  the  six  channels  with  respect  to  P    ,  that  is,  with 

respect  to  the  linear  J-scale.   Figure  10,  the  first  derivative 

with  respect  to  J  of  the  transmittances  in  Fig.  9,  is  the 
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Figure  9.   Carbon  dioxide  transmittances  channels  1  through  6, 
after  McMillin  et  al.  (1973. 


weighting  function  used  in  determining  this  weighted  mean 
Planck.   The  usual  procedure  as  used  by  Smith  et  al.  (19  72), 
involves  the  use  of  a  constant  normalized  wave  number 
(v  =  700  cm"  )  in  Eq.  37 ,  whereas  in  this  thesis,  a  normalized 
wave  number,  decreasing  with  decreasing  values  of  pressure  was 
found  to  give  a  more  accurate  result,  similar  to  results  of 
Moran  (1974) .   By  using  this  varying  wave  number  (Table  4) , 
advantage  can  be  taken  of  the  level  of  maximum  contribution  of 
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Table  4.   Normalized  wave  numbers  and  corresponding  pressure 
levels  used  in  inverting  the  Planck  function  for 
temperature  retrieval. 


Pressure 
mb 


P(J) 


Wave  Number  (v) 
cm""l 
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Figure  10.   Weighting  functions  used  in  determining  a 
weighted  mean  Planck  function. 
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that  channel  to  the  derived  equivalent  blackbody  temperature. 

The  weighted  mean  Planck  function  then  becomes 

6    (n) 

Z  B)n)   (J)  Ax. (J) 

,_.  i=l 

BW[v,T(J)]  =  .  (37) 

6 

Z   At. (J) 

i=l 


By  using  this  weighted  Planck,  the  100  mean  temperatures, 
T(J) ,  can  be  determined  by  the  inversion  of  Eq.  29 


T(J)  =  _f  V  (n)  _ •  OS) 

ClvJ  +  Bl  ;  (v,J) 

ln[- 


B(n)(v,J) 


In  Eqs .  37  and  38,  which  are  accomplished  after  the  iteration 
procedure  in  each  channel  is  complete,  the  values  of  At^(J) 
are  centered  differences  taken  over  A J  =  2 ,  but  centered  on 
each  individual  J  level.   A  forward  At.  difference  is  used  to 
initiate  the  process  at  level  J  =  1  and  a  backward  difference 
is  used  to  terminate  the  weighting  process  at  J  =  100.   In 
this  solution  procedure,  Eqs.  37,  38,  the  algorithm  used  is 
essentially  after  Smith  et  al.  (1972) ,  except  for  the  use  of 
the  variable  v . 

C.   STEPWISE  TEMPERATURE  RETRIEVAL 

The  procedure  for  retrieving  the  layer-mean  temperatures 
for  layer  J  =  1,2,..., 100,  from  the  clear-column  radiances 
can  be  summarized  as  follows  (program  names  refer  to  the 
listing  at  the  end  of  this  thesis) . 
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1.  Input  the  necessary  data. 

a.  sounding  latitude. 

b.  zenith  angle. 

c.  15-level  first  guess  temperatures. 

d.  sea  surface  temperature. 

e.  clear-column  radiances  in  channels  1  to  6 . 

2.  Expand  the  15-level  first  guess  temperature 
profile  to  56  levels  and  merge  it  with  climatology  at  the 
10  mb  level  using  subroutine  UPPER. 

3.  Expand  the  56-level  temperature  profile  to  100 
levels  using  Eq.  2  7  and  subroutine  LEVELL. 

4.  Calculate  the  atmospheric  transmittance  due  to 
carbon  dioxide  in  channels  1  through  6  at  each  level  J  using 
subroutine  RDTEMP  and  adjust  for  the  zenith  angle  by  using 
subroutine  MARTAU. 

5.  Calculate  the  water  vapor  content  of  the  atmosphere 
in  the  form  of  a  mixing  ratio  profile,  one  value  at  each  level 
J,  using  function  GOFF  and  Eqs.  7  and  8. 

6 .  Use  the  mixing  ratios  to  calculate  the  path 
corrected  atmospheric  transmittance  due  to  water  vapor  in 
channels  1  through  6  at  each  level  J  using  subroutine  TRANW. 
(Note:   steps  5  and  6  were  performed  on  the  first  iteration 
only,  due  to  the  small  change  in  water  vapor  transmissivity 
from  the  first  guess  to  the  final  temperature  profile. 
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7.  Calculate  the  total  atmospheric  transmissivity 
in  each  channel  and  at  each  level  through  the  product  of 
Eq.  3  which  uses  path  corrected  values  of  transmittance  due 
to  carbon  dioxide,  water  vapor,  and  ozone. 

8.  For  the  first  iteration,  using  the  100-level 
first  guess  temperatures,  the  Planck  functions  must  be 
calculated  in  each  of  the  6  channels. 

a.  Calculate  the  100  J-level  Planck  functions 
using  Eq.  29. 

b.  Calculate  the  33  K-level  Planck  functions  by 
fitting  a  fourth  order  Lagrangian  polynomial  to  the  100  J-level 
Planck-profile  [after  Conte  and  deBoor  (1965) ] . 

9.  Calculate  the  At (J)  for  each  channel  in  the  100 
J-level  and  At(K)  for  the  33  K-level  atmospheres.   The  At (K) 
will  be  used  in  the  discretized  form  (Eq.  30)  of  the  RTE ,  in 
order  to  reduce  the  truncation  error  associated  with  the 
finite  summation,  substituted  for  the  Riemann  integration  of 
the  integrand  function  in  Eq.  28.   This  is  important  while 
dealing  with  the  t(J)  values  which  occur  before  the  temperature 
profile  and  its  associated  transmissivities  have  converged. 
The  converged  At (J)  will  be  used  with  Eq.  37  in  producing  the 
new  100-level  temperature  profile. 

10.  Compute  the  layer-mean  Planck  values  by  Eq.  31. 

11.  Compute  the  radiance  in  each  channel  using  Eq.  30 
and  determine  the  residual  difference  between  the  clear-column 
radiance  and  the  computed  radiance,  Eq.  33. 
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12.  If  the  residual  value  is  not  close  to  zero, 
another  iteration  must  be  performed  to  improve  the  accuracy 
of  the  retrieved  temperature  profile. 

a.  Calculate  the  now  100-level  and  33-level 
Planck  functions  by  adding  the  residual  to  each  previous 
value,  Eq .  34. 

b.  Calculate  the  new  100-level  temperature 
profile  using  Eqs.  37  and  38,  and  go  back  to  step  4.   Continue 
recycling  until  convergence  is  achieved  (see  Eq.  36) . 

13.  If  the  residual  is  small  enough  (Eq.  36  is  true), 
go  directly  to  the  calculation  of  the  new  100-level  temperature 
profile,  Eqs.  37  and  38. 

D.   PROGRAM  OUTPUT 

The  program  is  designed  to  produce  a  mean  layer  temperature 
profile  at  100  levels  from  1000  mb  to  0.01  mb.   For  the  purpose 
of  testing  the  accuracy  of  these  retrieved  temperatures, 
additional  inputs  were  needed. 

First  a  set  of  radiosonde  soundings,  as  close  as  possible 
to  each  retrieval  site  in  latitude,  longitude,  and  real  time 
was  read  in.   While  these  soundings  are  not  considered  to  be 
the  exact  state  of  the  atmosphere  at  the  retrieval  site,  they 
are  as  close  as  possible  and  therefore  were  used  to  determine 
the  error  in  the  retrieval  program  at  each  level  J  by  using 

ERROR(J)  =  T(J)  -  AT  (J)  (39) 
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where   T(J)  is  the  retrieved  temperature,  and 

AT (J)  is  the  radiosonde  sounding  temperature. 

Secondly,  a  set  of  retrieved  temperatures  supplied  from 
NESS  for  each  case  investigated  was  read  in  and  compared  with 
the  radiosonde  temperatures . 

The  program  then  output  the  following  for  each  of  the 
2  33  retrievals  attempted  (see  sample  output  following  the 
program  listing) : 

1.  Station  number  (1,2, . . . ,233) 

2.  Sounding  zenith  angle  in  degrees. 

3.  Station  latitude. 

4.  The  adjustment  at  10  mb  necessary  to  fit  the  first 
guess  temperature  profile  to  the  climatological 
temperatures . 

5.  The  number  of  iterations  necessary  for  convergence. 

6.  Table  of  results,  listing  for  each  of  the  15  levels: 

a.  Pressure 

b.  Retrieved  temperature 

c.  Radiosonde  sounding  temperature 

d.  Error  (retrieved  --  sounding) 

e.  NESS  retrieved  temperature 

f.  Error  (NESS  retrieved  -  sounding) 

In  order  to  easily  compare  the  accuracy  of  the  retrieval 
with  the  radiosonde  data  and  also  to  compare  the  error  of 
this  method  with  that  of  NESS,  histograms  were  drawn  and 
statistical  data  were  obtained.   The  errors  were  compared 
level  for  level  on  an  overall  basis  (globally  averaged)  and 
also  on  a  latitude-band  basis  for  the  latitude  ranges  shown 
in  Table  5. 
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Table  5.   Latitude  bands  used  for  error  analysis  of  the 
VTPR  retrieved  temperatures. 


LATITUDE  BAND  BAND  WIDTH 

1  0  <  LAT  <  20N 

2  20N  <_  LAT  <  4 ON 

3  40N  <_  LAT  <  60N 

4  60N  <  LAT 
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V.   RESULTS  OF  THE  VTPR  TEMPERATURE  RETRIEVAL 

In  general,  the  errors  (retrieved  temperature  -  radiosonde 
temperature)  were  less  than  5  °K  when  averaged  over  the 
latitude-bands  shown  in  Table  5.   The  errors  also  compared 
favorably  with  the  results  of  the  NESS  VTPR  retrieval  in  all 
15  atmospheric  levels  and  the  four  latitude  bands.   In 
calculating  the  errors,  the  assumption  that  the  radiosondes 
matched  the  VTPR  retrieval  sites  exactly  in  time  and  space 
was  made.   The  difference,  though  small  (averaging  ±0.6 
degrees  in  latitude,  ±0.7  degrees  in  longitude,  and  ±1.87 
hours  in  time) ,  coupled  with  the  small  data  base  of  2  33 
retrievals,  could  have  had  some  effect  on  the  results  and 
should  be  of  some  consideration. 

A.   GLOBAL  RESULTS 

The  global  mean  results  obtained  here  are  summarized  in 
Table  6  along  with  the  averaged  results  from  the  NESS 
retrievals.   The  frequency  distributions  of  the  global  error 
analysis  at  each  of  15  key  atmospheric  levels  are  shown  as  a 
function  of  the  error  classes  encountered  (See  Appendix  A-l 
to  A-15) .   These  frequency  distributions,  in  the  form  of 
histograms,  are  included  as  a  visual  aid  in  analyzing  the 
deviation  of  the  error  from  the  mean.   They  show  that  even 
though  the  range  of  errors  can  extend  over  as  much  as  20  °K, 
the  probability  of  a  relatively  large  error  is  much  less  than 
the  probability  of  the  small  errors  in  the  vicinity  of  the  mean 
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Table  6.   Retrieved  temperature  error  analysis  on  a  global  basis 


RETRIEVAL 

NE! 

3S 

PRESSURE    (mb) 

MEAN 

RMS  ERROR 

MEAN 

RMS  ERROR 

NO.    OF 

ERROR    (°K) 

ERROR    (°K) 

SAMPLES 

10 

-0.55 

4.00 

-0.75 

4.51 

59 

20 

0.57 

3.25 

-0.11 

2.72 

112 

30 

-0.06 

3.20 

-0.17 

3.12 

142 

50 

0.37 

3.20 

-0.17 

3.12 

142 

70 

0.75 

4.39 

0.91 

5.12 

179 

100 

-0.45 

3.89 

-2.30 

5.12 

179 

150 

-0.30 

2.58 

-1.51 

2.73 

208 

200 

-0.66 

3.33 

-1.04 

3.29 

215 

250 

0.15 

3.22 

0.41 

3.33 

206 

300 

-0.58 

2.82 

0.24 

2.70 

223 

400 

-0.18 

2.79 

-0.09 

2.07 

226 

500 

-0.67 

3.09 

-0.32 

2.13 

227 

700 

-0.99 

3.29 

-0.86 

2.81 

226 

850 

-0.08 

3.26 

-0.31 

3.55 

227 

1000 

-0.24 

1.77 

0.26 

3.89 

196 
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B.   LATITUDE  BAND  RESULTS 

Because  of  space  limitations,  frequency  distributions 
were  not  included  for  the  latitude  band  results.   However, 
Table  7  furnishes  a  summary  of  the  latitude  band  analysis 
of  errors  including  their  root  mean  square  errors  along  with 
the  corresponding  NESS  results. 
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VI.   CONCLUSIONS 

The  first  objective  of  this  thesis  was  to  incorporate  the 
state-of-the-art  methods  [as  used  by  Fritz  et  al .  (1972)]  in 
calculating  atmospheric  transmittances ,  into  the  retrieval 
scheme.   The  method  involved  here  for  calculating  the  effects 
of  water  vapor,  ozone,  and  carbon  dioxide  transmittance  proved 
to  be  satisfactory.   This,  coupled  with  the  inclusion  of  a 
33  level  integration  scheme  and  the  use  of  a  variable  wave 
number,  as  opposed  to  the  constant  v"  =  700  cm   ,  combined  to 
give  results  comparable  to  those  of  NESS. 

The  temperature  profile  convergence  occurred  after  an 
average  of  5  iterations.   While  this  is  an  improvement  over 
the  Moran  (1974)  model,  as  would  be  expected  with  the  more 
sophisticated  procedures  incorporated,  the  number  of  iterations 
required  would  possibly  be  reduced  further  through  the  use  of 
a  more  accurate  first  guess  temperature  profile.   Gelman  et 
al.  (1972)  used  a  regression  technique  for  estimating  the 
first  guess  temperatures  above  10  mb,  as  opposed  to  merging 
with  the  standard  atmosphere  climatology  as  used  here.   Since 
this  regression  technique  of  Gelman  produces  a  more  accurate 
first  guess  temperature  profile,  it  would  reduce  the  number 
of  iterations  required  for  convergence  even  further.   However 
this  program  was  not  available  locally  for  inclusion  in  this 
thesis  . 
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The  second  objective,  that  of  evaluation  through  compari- 
son with  matching  radiosondes,  showed  that  the  retrieval 
system  produced  errors  that  were  comparable  to  NESS  retrievals, 
or  better,  especially  in  the  area  of  the  tropopause.   The 
small  data  base  of  233  soundings  is  not  considered  large 
enough  to  evaluate  the  accuracy  of  the  retrieval  system  beyond 
that,  and  more  experimentation  with  a  larger  sample,  and 
inclusion  of  the  Gelman  regression  technique  is  required. 

It  is  concluded  that  this  thesis  has  demonstrated  a 
consistent  and  accurate  means  of  temperature  profile  retrieval 
from  VTPR  clear-column  radiance  data. 
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APPENDIX  A-l 


Histogram;of  the  globally  averaged  errors  showing  probability 
distribution  [P(E)]  of  the  occurrence  of  the  errors,  the  actual 

^hr°M     ?m? n°ie^  ^  thG  err°r  Class  E  <marked  along  the  top)  , 
the  Mean  (M) ,  and  the  Quartiles  (.),  for  sample  size  196 
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APPENDIX  A-2 


Histogram  of  the  globally  averaged  errors  showing  probability 
distribution  [P(E)]  of  the  occurrence  of  the  errors,  the  actual 
error  frequencies  in  the  error  class  E  (marked  along  the  top) , 
the  Mean  (M) ,  and  the  Quartiles  (•)/  for  sample  size  227. 
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APPENDIX  A-3 

Histogram  of  the  globally  averaged  errors  showing  probability 
distribution  [P(E)J  of  the  occurrence  of  the  errors,  the  actual 
error  frequencies  in  the  error  class  E  (marked  along  the  top) , 
the  Mean  (M) ,  and  the  Quartiles  (•),  for  sample  size  226. 
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APPENDIX  A-4 

Histogram  of  the  globally  averaged  errors  showing  probability 
distribution  [P(E)]  of  the  occurrence  of  the  errors,  the  actual 
error  frequencies  in  the  error  class  E  (marked  along  the  top) , 
the  Mean  (M) ,  and  the  Quartiles  (•) /  for  sample  size  227. 
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APPENDIX  A- 5 

Histogram  of  the  globally  averaged  errors  showing  probability 
distribution  [P(E)]  of  the  occurrence  of  the  errors,  the  actual 
error  frequencies  in  the  error  class  E  (marked  along  the  top) , 
the  Mean  (M) ,  and  the  Quartiles  (•),  for  sample  size  226. 
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APPENDIX  A-6 

Histogram  of  the  globally  averaged  errors  showing  probability 
distribution  [P(E)]  of  the  occurrence  of  the  errors,  the  actual 
error  frequencies  in  the  error  class  E  (marked  along  the  top) , 
the  Mean  (M) ,  and  the  Quartiles  (•) ,  for  sample  size  223. 


FREQUENCIES 
3    1    2 


20       31       42 


35      37 
—  + +  ■ 


20       14 
+ *- 


P(E) 

.18 

.16 

.14 
.12 

.10 
.08 
.06 

.04 
.02 


*** 


***+**  *4-«  ##+•**: 


-Jt-  <JL- 

ju  ju 


* 

V  3 
*: 


* 
* 
* 
* 
* 
+  * 


-9.  19 


-6.82 


-4.45 


+ 

I 

-2. 


*#*, 

,M 

*** 

,M 

#*#, 

»M 

#♦#, 

,M 

###  , 

,M 

*£*, 

►  M 

***, 

,M 

V 

.* 

***, 

,M 

.    ^ 

*r  *i*  *r  , 

,  M** 

.* 

#*#, 

,  M#* 

,* 

**#, 

M** 

* . 

,    * 

*  *  * , 

,  N|£# 

* 

,* 

***, 

,  M## 

'i* 

,    * 

^•,  o,  *v 

,  M** 

* 

.    * 

#£#, 

,M** 

t    ^ 

##*, 

,  M** 

V 

,    * 

***, 

,  M5^"^ 

>* 

***, 

,  M** 

1*   , 

,* 

##*, 

,  M** 

%'• 

►  * 

#*#  , 

,  M** 

,    * 

•*c-  -c-  *c-  t 

,M*« 

s'- 

1* 

*##  . 

M*# 

=:c 

.    * 

i-  'i*  i*  , 

***■  i 

t         '.- 

*  #*, 

,  M** 

;'; 

,* 

*#*   , 

,  M*# 

*! 

>* 

•A.   ^  -X, 

M** 

*r  ( 

,* 

♦  ##, 

*  < 

,     * 

£**, 

,  M#* 

JL. 

,  V 

#**    , 

.  M## 

•IT    , 

.  * 

#**< 

.  M** 

^    1 

VI-     JU     V* 
*T     *V*     ***       ( 

,M** 

*| 

.* 

.  M** 

♦ 

,  * 

•v- 1*  "V"  , 

,  M*# 

v(- 

£#* , 

,  M** 

'O    ( 

!  * 

***, 

,M** 

*, 

,  * 

&##  , 

.  M** 

*i 

* 

#  'r^  4 

,  M  *  * 

* 

**#, 

»  M*  * 

Jo- 

,  * 

¥    V    1*    , 

,  M** 

JU 

1*   | 

* 

JL-     -JV  «/~ 

•^r  tt  V  4 

,  M*=* 

*■ 

£ 

*## , 

,M** 

*  . 

,  * 

#  #&  < 

.M*# 

*. 

.* 

**#  , 

,  M** 

*  , 

,  V 

***, 

M#* 

*. 

* 

#*#, 

M** 

mC 

* 

#**, 

M*# 

.* 

>r  =£  *  4 

M** 

* 

ji  ^c^t 

,M*« 

*    . 

.  * 

**#] 

M*# 

#« 

* 

***, 

,  m**  +  * , 

08 


vt,    -J-  JL. 

JL-    -JU  JU 


0.29 


>'f     JL. 

I 
2. 


5?  -r 

* 


3p  5p  *p. 

*** 
*** 
*** 

•.LO.  Jf 

-,»    -y.     ,,i 

;,x  ?J<  -V 


* 


^-  V  5JC 


;=* 


66 


5.03 


-1?  4-  ^r  ■**  * 

7.40 


ERROR  CLASSES  <E) 
300  MB   HISTOGRAM  FOR  GLOBALLY  AVERAGED  ERRORS 


71 


APPENDIX  A-7 

Histogram  of  the  globally  averaged  errors  showing  probability 
distribution  [P(E)]  of  the  occurrence  of  the  errors,  the  actual 
error  frequencies  in  the  error  class  E  (marked  along  the  top) , 
the  Mean  (M) ,  and  the  Quartiles  (•),  for  sample  size  206. 
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APPENDIX  A- 8 

Histogram, of  the  globally  averaged  errors  showing  probability 
distribution  [P(E)]  of  the  occurrence  of  the  errors,  the  actual 
error  frequencies  in  the  error  class  E  (marked  along  the  top) , 
the  Mean  (M) ,  and  the  Quartiles  (•),  for  sample  size  215. 
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APPENDIX  A- 9 

Histogram, of  the  globally  averaged  errors  showing  probability 
distribution  [P (E) ]  of  the  occurrence  of  the  errors,  the  actual 
error  frequencies  in  the  error  class  E  (marked  along  the  top) , 
the  Mean  (M) ,  and  the  Quartiles  (•)/  for  sample  size  208. 
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APPENDIX  A-10 

Histogram  of  the  globally  averaged  errors  showing  probability 
distribution  [P (E) ]  of  the  occurrence  of  the  errors,  the  actual 
error  frequencies  in  the  error  class  E  (marked  along  the  top) , 
the  Mean  (M) ,  and  the  Quartiles  (•)/  for  sample  size  179. 
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APPENDIX  A-ll 

Histogram  of  the  globally  averaged  errors  showing  probability 
distribution  [P(E)]  of  the  occurrence  of  the  errors,  the  actual 
error  frequencies  in  the  error  class  E  (marked  along  the  top) , 
the  Mean  (M) ,  and  the  Quartiles  (•),  for  sample  size  179. 
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APPENDIX  A-12 

Histogram  of  the  globally  averaged  errors  showing  probability 
distribution  [P(E)]  of  the  occurrence  of  the  errors,  the  actu 
error  frequencies  in  the  error  class  E  (marked  along  the  top) 
the  Mean  (M) ,  and  the  Quartiles  (•),  for  sample  size  142 
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APPENDIX  A-13 

Histogram  of  the  globally  averaged  errors  showing  probability 
distribution  [P(E)   of  the  occurrence  of  the  errors,  the  actual 

XrS     ?MfnC1^  ln  thG  err°r  ClaSS  E  (marked  al°^  the  top?, 
the  Mean  (M) ,  and  the  Quartiles  (•),  for  sample  size  142 
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APPENDIX  A-14 


Histogramof  the  globally  averaged  errors  showing  probability 
distribution  [P(E)]  of  the  occurrence  of  the  errors,  the  actual 
error  frequencies  in  the  error  class  E  (marked  along  the  top) , 
the  Mean  (M)  ,  and  the  Quartiles  (•),  for  sample  size  112. 
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APPENDIX  A-15 

Histogram  of  the  globally  averaged  errors  showing  probability 
distribution  [P(EJ]  of  the  occurrence  of  the  errors,  the  actual 
error  frequencies  in  the  error  class  E  (marked  along  the  top) , 
the  Mean  (M) ,  and  the  Quartiles  (•),  for  sample  size  59 
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APPENDIX  B 


OZONE  TRANSMITTANCES  FOP  THE  SIX  VT2*  CLEAR-COLUMN  CHANNELS 
LEVEL   CHAN  1     CHAN  ?  CHAN  3     CHAN  4     CHAM  5     CHAN  6 


1  0.999992 

2  0.999982 

3  0.999966 

4  0.999941 

5  0.999906 

6  0.999859 

7  0.999797 

8  0.999717 

9  0.999618 

10  0.999497 

11  0.999352 

12  0.999182 

13  0.998982 

14  0.998748 

15  0.998476 

16  0.998162 

17  0.997301 

18  0.997392 

19  0.996937 

20  0.9964-40 

21  0.995905 

22  0.995331 

23  0.994703 

24  0.993998 

25  0.993196 

26  0.992294 

27  0.991320 

28  0.990306 

29  0.989287 

30  0.988294 

31  0.987361 

32  0.986509 

33  0.985736 

34  0.985037 

35  0.984405 

36  0.983835 

37  0.983323 

38  0.982861 

39  0.982445 

40  0.982069 


0.999987 
0.999974 
0.999953 
0.999917 
0.999865 
0.999795 
0.999703 
999586 
999440 
999263 
999051 
998602 
0.993510 
0.998172 
0.997781 
0.997333 
0.996823 
0*996249 
0.995611 
0.994915 
0.994163 
0.993353 
0.992469 
0.991490 
0.990395 
0.989178 
0.987867 
0.986489 
0.985076 
0.983657 
0.982262 
0.980911 
0.979606 
0.978345 
0.977129 
0.975954 
0.974821 
0.973728 
0.9726  74 
0.971658 


0.999986 
0.999972 
0.999949 
0.999914 
0.999366 
0.999804 
0.999725 
0.999626 
0.999504 
0.999358 
0.9991 83 
0.998979 
0.998742 
0.998466 
0.9981 48 
0.997783 
0.997365 
0.996392 
0.996361 
0.99  5  7  73 
0.995124 
0.994415 
0.993645 
Oc9923  12 
0.991916 
0.990952 
0.939911 
0.988731 
0.937552 
0.986212 
0.984751 
0.983165 
0.981446 
0.979671 
0.977793 
0.975349 
0.973853 
0.9718  21 
0.969768 
0.967703 


0. 999987 
0.99  9  96  3 
0.999940 
3.999900 
0.999  8  44 
0.999771 
0.99  96  79 
0.999565 
0.999425 
0.999253 
0.999059 
0.993326 
0.9  93  5  54 
0.998239 
0.997374 
0.997453 
0.996973 
0.995429 
0.995823 
0.995156 
0.9  9  44  30 
0.99  3  643 
0.99  2  7  83 
0.991334 
0.993779 
0.989614 
0.93R352 
0.937012 
0.935610 
0.934165 
0.98  2  693 
0.981 209 
0.979717 
0.973222 
0.976728 
0.975240 
0.973761 
0.972295 
0.970843 
0.969423 


0.999991 
0.999974 
0.999949 
0.999^12 
0.999857 
0.999780 
0.999679 
0.999551 
0.999392 
0.999198 
0.998967 
0.998694 
0.99  83  75 
0.998306 
0.997579 
0.997092 
0.996537 
0.995914 
0.995224 
0.9944  70 
0.993656 
0.992781 
0.991827 
0.990771 
0.989591 
0.938283 
0.986372 
0.985391 
0.983869 
0.932333 
0.983829 
0.979362 
0.977941 
0.976565 
0.975231 
0.973939 
0.972689 
0.971479 
0.970309 
0.969177 


0.999992 
0.999935 
0.999971 
0.999943 
0.999914 
0.999369 
0.999813 
0.99  97  3  4 
0.999643 
0.9995  2  5 
0.999386 
0.999223 
0.999032 
0.998313 
0.993553 
0.993253 
997922 
997543 
997121 
996^7 
996152 
995506 
995012 
0.994364 
0.993655 
0.9  9288? 
0.992053 
0.991162 
0.990224 
0.939240 
0.988214 
0.987153 
0.98506? 
0.984945 
0.983310 
0.982662 
0.981507 
0.980350 
0.^73197 
0.978054 
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0.947737 
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99 

0.972478 
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0.947564 
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00 

0.972376 

0.952297 

0.939628 

0.945335 

0.947391 

0.960953 
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CARBON  DIOXIDE  TR ANSMI TTANCES  FO*  THE  SIX  VTP*  CHANNELS 


VEL 

CHAN    1 

CHAN    2 

CHAN    3 

r 
I 

;ha\i   4 

r 
I 

;han  5 

CHANJ    6 

1 

0. 

991131 

0, 

.999653 

0,999886 

0. 

9^991 ? 

0. 

999936 

0, 

,999993 

2 

0. 

988616 

0, 

,999418 

0.999733 

0. 

999799 

0. 

999842 

0. 

,999935 

3 

0. 

985314 

0, 

.999092 

0.9^9537 

0. 

999575 

0. 

999653 

0, 

,999962 

4 

0. 

980754 

0, 

.9986  08 

0.999259 

0. 

999173 

0. 

999  3  08 

0. 

,999939 

5 

0. 

975299 

0. 

.997973 

0.998884 

0. 

998641 

0. 

998848 

0. 

,999315 

6 

0. 

969083 

0, 

.997204 

0.998390 

0. 

998053 

0. 

,998348 

0, 

,99967* 

7 

0. 

96  1684 

0, 

.996247 

0.997750 

0. 

997415 

0. 

997822 

0. 

,999*85 

8 

0. 

952678 

0. 

.995045 

0.996952 

0, 

,996  73  8 

0. 

,Q97281 

0. 

,999253 

9 

0. 

941640 

0. 

.993540 

0.995971 

0. 

9  9  6  D  2  2 

0. 

,996731 

0. 

,993991 

10 

0. 

928429 

0, 

.991679 

0.994789 

0. 

995254 

0. 

9961 55 

0. 

,998697 

11 

0. 

913059 

0, 

,989418 

0.993404 

0. 

9  9  *  4  3  3 

0. 

995537 

0. 

,998403 

12 

0. 

895530 

0, 

,986719 

0.991792 

0. 

993556 

0. 

,994859 

0. 

,  998  101 

13 

0. 

875965 

0, 

.983537 

0.989937 

0. 

,992624 

0. 

,99  41 12 

0. 

,907802 

14 

0. 

854614 

0, 

.979809 

0.987313 

0. 

99162  7 

0. 

,993282 

0. 

,997503 

15 

0. 

831837 

0. 

,975464 

0.985379 

0. 

990533 

0. 

992342 

0. 

997193 

16 

0. 

808120 

0. 

,970429 

0.982596 

0. 

9  B  9  3  2  4 

0. 

,991280 

0. 

,996865 

17 

0. 

783978 

0. 

,964630 

0.979*?6 

0. 

987985 

0. 

990083 

0. 

,996519 

18 

0. 

7598  06 

0, 

.957937 

0.975819 

0. 

9  36  496 

0. 

988739 

0, 

,996145 

19 

0. 

735895 

0, 

,950424 

0.971 736 

0. 

98*8  47 

0. 

987244 

0. 

.995733 

20 

0. 

712363 

0, 

,941861 

0.967135 

3. 

983030 

0. 

,935595 

o, 

9953  01 

21 

0. 

689153 

0. 

,932223 

0.961984 

0. 

981041 

0. 

083792 

5, 

,99*323 

22 

0. 

666067 

0, 

.921432 

0.956255 

0. 

973881 

0. 

,981841 

0, 

,994319 

23 

0. 

642342 

0, 

,909414 

0.949899 

0. 

9  7  6  5  31. 

0. 

979731 

0. 

,993767 

24 

0. 

619215 

0, 

,896098 

0.942897 

0, 

,973993 

0. 

,977470 

0. 

,993173 

25 

0. 

594999 

0, 

,831417 

0.9352 16 

0. 

9  712  5  3 

0. 

,975059 

0, 

,9925*3 

26 

0. 

570067 

0, 

.865308 

0.926777 

0. 

968284 

0. 

972*73 

0, 

.991863 

27 

0. 

544428 

0. 

,847723 

0.91 7538 

0. 

9S5057 

0. 

,969713 

0. 

,991121 

28 

0. 

518133 

0, 

.828618 

0.907477 

0. 

961 575 

0. 

,966787 

0. 

,990321 

29 

0. 

491282 

0, 

.807966 

0.896576 

0. 

957833 

0. 

,953705 

0, 

.989*65 

30 

0. 

463997 

0. 

.785753 

0.884823 

0. 

,953848 

0. 

,960479 

0, 

,988555 

31 

0. 

436414 

0, 

,761986 

0.872210 

0. 

949607 

0. 

957125 

0. 

,987592 

32 

0. 

408677 

0. 

,73o687 

0.858738 

0. 

,945123 

0. 

,953658 

0. 

,936582 

33 

0. 

380939 

0, 

.709907 

0.844414 

0. 

,9*0  40  8 

0. 

,950097 

0. 

,985527 

34 

0. 

353356 

0, 

.681707 

0.829254 

0. 

935474 

0. 

9*6*59 

0. 

,984435 

35 

0. 

326084 

0, 

.652176 

0.813267 

0. 

,930328 

0. 

,942757 

0. 

,983338 

36 

0. 

299286 

0, 

.621431 

0.796393 

0. 

,924929 

0. 

,933977 

0. 

,982135 

37 

0. 

273120 

0. 

,589637 

0.778590 

0. 

,919237 

0, 

,935107 

0. 

.98090* 

38 

0. 

247735 

0, 

.556957 

0.759887 

0. 

,913260 

0. 

,9311 57 

0, 

,979613 

39 

0. 

223274 

0. 

.523557 

0.740321 

0. 

907008 

0. 

,927138 

0. 

,97328? 

40 

0. 

199827 

0, 

.489605 

0.720067 

0. 

,900580 

0. 

,923095 

0, 

,976925 
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41 

42 
43 
44 
45 
46 
47 
48 
49 
50 
51 
52 
53 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
66 
67 
68 
69 
70 
71 
72 
73 
74 
75 
76 
77 
78 
79 
80 


0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 


177535 

156501 

136779 

118460 

101568 

086144 

072217 

059791 

048830 

039287 

031120 

024258 

018534 

0139  70 

010292 

007424 

005227 

003538 

002388 

CO  1541 

000968 

000591 

000341 

000180 

000082 

000020 

000000 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0.4 
0.4 
0.3 
0.3 
0.3 
0.2 
0.2 
0.2 
0.1 
0.1 
0.1 
0.1 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 


55338 
20970 
86711 
52818 
19468 
86895 
55350 
25062 
96242 
69099 
43870 
20747 
99830 
81190 
64897 
50927 
39151 
29433 
216  29 


0.010836 

0.007353 

0.0048  10 

0.003047 

0.001836 

0.001156 

0.000691 

0.000395 

0.000221 

0.000123 

0.000060 

0.000019 

0.000000 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


0.699027 
0.677231 
0.654787 
0.6316  11 
0.607929 
0.583764 
0.559096 
0.533982 
0.508515 
0.482338 
0.456970 
0.430848 
0.404528 
0.373156 
0.351884 
0.325359 
0.3002  36 
0.275213 
0.250973 

ri  -i  ->  -7  c  i  "» 
U  .    C  C   I    O    L  C 

0.205211 
0.133853 
0. 163613 
0.144515 
0.126588 
0.109356 
0.094406 
0.380334 
0.06  76  55 
0.056346 
0.046403 
0.037766 
0.030335 
0.023984 
0.018661 
0.014292 
0.010759 
0.007967 
0.005844 
0.004220 


0.89 


88 
87 
87 
86 
85 


0.84 
0.83 
0.8? 
0.82 
0.81 
0.79 
0.7Q 
0.77 
0.76 
0.75 
0.73 
0.72 
0.70 
0.60 
0.6  7 
0.6  5 
0.53 
0.61 
0.5  9 
0.57 
0.55 
0.52 
0.53 
0.47 
0.  45 
0.42 
0.  39 
0.37 
0.34 
0.31 
0.29 
0.2S 
0.23 
0.21 


3880 
6900 
9668 
?086 

4  2  6  8 
6180 
7746 
8040 

9  7  70 
0248 
0294 
9742 
3  5  02 
6  5  50 
3863 
042  7 
6  2  3  5 
1332 

5  756 

2462 
4664 
S031 
6461 
5843 
4072 
120  2 
7442 
?S86 
7534 
1509 
5025 
3223 

10  34 
38  2  5 
6940 
0388 
4438 
9271 
5037 


0.91 
0.91 
0.91 
0.90 

0.90 
0.89 
0.89 
0.88 
0.83 
0.87 
0.87 
0.85 
0.86 
0.85 
0.8  4 
0.84 
0.83 
0.82 
0.82 

On  « 
.  O  t 

0.80 
0.79 
0.78 
0.77 
0.75 
0.75 
0.74 
0.73 
0.72 
0.70 
0.69 
0.68 
0.66 
0.64 
0.63 
0.61 
0.59 
0.58 
0.56 
0.54 


8934 
4800 
0554 
6201 
1789 
7299 
2698 
7974 
31  25 
3152 
3021 
7663 
2036 
6124 
9  9  09 
3377 
65  17 
9334 
1833 

-» :*  t>  c 

5783 
7179 
8142 
8612 
8520 
7793 
6433 
4486 
1936 
3741 
4903 
0453 
5387 
9593 
3165 
6200 
3623 
0505 
18  78 
2755 


0.975523 
0.974065 
0.972570 
0.971004 
0.969  40? 
0.967754 
0.966037 
0.964244 
0.962376 
0.960436 
0.958  401 
0.956221 
0.953867 
0.951329 
0.948595 
0.945657 
0.942505 
0.939145 
0.93558? 
0.93179  6 
0.927765 
0.923464 
0.913861 
0.913897 
0.903503 
0.902609 
0.896197 
0.889290 
0.831867 
0.8738  71 
0.855283 
0.856134 
0.846389 
0.835938 
0.824852 
0.813206 
0.800919 
0.783043 
0.774598 
0.760583 


110 


81 

0.0 

0.0 

0.002995 

0.  191390 

0.523149 

0.745  9  95 

82 

0.0 

0.0 

0. 002071 

0.  170036 

0.503112 

0.73087'+ 

83 

0.  0 

0.0 

0.001380 

0.14965^ 

0.482760 

0.71530? 

84 

0.0 

0.0 

0.030385 

0. 130787 

0.452181 

0.699323 

85 

0.0 

0.0 

0.000548 

0.113420 

0.441413 

0.68295+ 

86 

0.0 

0.0 

0.000334 

0.09  7  54  7 

0.420495 

0.666182 

87 

0.0 

0.0 

0.000208 

0.083170 

0.399466 

0.649015 

88 

0.0 

0.0 

0.000132 

0.070295 

0.378367 

0.631455 

89 

0.0 

0.0 

0.0 

0.058915 

0.357233 

0.61*499 

90 

0.0 

0.0 

0.0 

0.048942 

0.336122 

0.595161 

91 

0.0 

0.0 

0.0 

0.040301 

0.315132 

0.576553 

92 

0.0 

0.0 

0.0 

0.032900 

0.294559 

0.557781 

93 

0.0 

0.0 

0.0 

0.0?5612 

0.274311 

0.538379 

94 

0.0 

0.0 

0.0 

0.021319 

0.254492 

0.519862 

95 

0.0 

0.0 

0.0 

0.015910 

0.235156 

0.500744 

96 

0.0 

0.0 

0.0 

0.013286 

0.216375 

0.481554 

97 

0.0 

0.0 

0.0 

0.010351 

0.198231 

0.462330 

98 

0.0 

0.0 

0.0 

0.008008 

0.180793 

0.443100 

99 

0.0 

0.0 

0.0 

0.005171 

0. 164126 

0.423879 

00 

0.0 

0.0 

0.0 

0.004732 

0. 148023 

0.404273 
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WATER  VAPOR  TRANSM ITT ANCES  FOR  THE  SIX  VTPR  CHANNELS 


LEVEL   CHAN  1 


CHAN  2 


CHAN  3 


CHAM  4 


CHAN  5 


CHAM  6 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 


1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 
0. 
0. 
0. 
0. 


000000 
000000 
000000 
000000 
000000 
OOCOOO 
000000 
000000 
000000 
000000 
000000 
000000 
000000 
OOOOOC 
000000 
000000 
000000 

cooooo 

000000 
000000 
000000 
000000 
000000 
000000 
000000 
000000 
000000 
000000 
000000 
000030 
000000 
000000 
000000 
000000 
000000 
000000 

999999 
999999 

999998 
999997 


000000 
000000 
000000 
000000 
000000 
000000 


1 .000000 


000000 
000000 
000000 
000000 
000000 
000000 
000000 
000000 
000000 
000000 
000000 
000000 
0  00  0  00 
000000 
000000 
000000 
000000 
000000 
000000 
000000 
000000 
000000 
000000 
000000 
000000 
000000 
000000 
000000 
000000 
000000 
000000 
1 .000000 
1.000000 


000000 

000300 
000000 
000000 
000000 
000000 
000000 
000000 
003000 
1.000003 
1. 0000  00 
1.000000 
1.000000 
1.000000 
1.000000 
1.000000 
1.000000 
1.000300 
1.000900 


000000 
000000 
000000 
000000 
0000  00 
000000 
000000 
000000 
000000 
OCOOOO 
000000 
000900 
000000 
000900 


1.000900 


000000 
009000 
000000 
0000  00 
000000 
000000 


1.09 
1  .09 
1.09 
1.09 
1.09 
1.09 
1.09 
1.03 
1.03 
1.09 
1.00 
1  .09 
1.00 
1.09 
1.00 
1.09 
1.09 
1  .09 
1.00 
1.09 
1  .09 
1.09 
1  .00 
1.03 
1.09 
1.09 
1.00 
1.03 
1  .03 
1.00 
1.03 
1  .03 
1.00 
1.03 
1.00 
1.03 
1  .03 
1.03 
1.00 
1.03 


0000 
3300 
9000 
9900 
9090 
9900 
9000 
9000 
9900 
0000 
0000 
9990 
0000 
0000 
0000 
9900 
9000 
9000 
9000 
9000 
9000 
9900 
0000 
9900 
9000 
0000 
0000 
3000 
0000 
0000 
0000 
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